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Abstract 

Relation of ultrametric analysis, wavelet theory and cascade models of turbulence is 
discussed. We construct the explicit solutions for the nonlinear ultrametric integral equation 
with quadratic nonlinearity. These solutions are built by means of the recurrent hierarchical 
procedure which is analogous to the procedure used for the cascade models of turbulence. 



1 Introduction 

p-Adic analysis and mathematical physics attract a lot of attention, see pQ, [2J, [3J. p-Adic 
pseudo differential operators (in particular, the Vladimirov operator D a ) and pseudodifferential 
equations (analogues of the equations of mathematical physics) were studied. p-Adic wavelets 
were introduced and the relation of wavelet analysis to spectral theory of p-adic pseudodifferential 
operators was described in [?]. Analysis of pseudodifferential operators and wavelets on general 
(locally compact) ultrametric spaces was developed in [5], [B], [7|. 

In the present paper we construct families of explicit solutions for the following nonlinear 
integro-differential equation: 

jU(x,t) + J J F(sap(x,Z,7i))v(£,t) (y(7i,t)-v(x,t))du(£)du(ri)+ 

+ J G(sup(x,Z)){v(x,t)-v((,t))dv(Z) = 0. 

Here the argument x G X belongs to the ultrametric space X, and t (time) is a real parameter. 
For the investigation of this equation we use the methods of ultrametric wavelet analysis, see the 
Appendices 1 and 3 for the discussion. 

We show that the above equation possesses solutions in the form of the product of ultrametric 
wavelet and exponent of time 

v(x,t) = e^ty/ifa)- 

Moreover, we construct also more general solutions for the equation above. We construct the 
unique solution of Cauchy problem with the initial condition in the space Dq(X) of mean zero test 
functions. These solutions are built by means of the recurrent hierarchical procedure. Procedures 



of this kind are related to the cascade models of turbulence. Using this observation we will call 
the equation above the quadratic cascade equation. Ultrametric (p-adic) analysis was applied for 
the first time to description of cascade models of turbulence in [8]. 
The structure of the paper is as follows. 

In Section 2 we construct solutions of Cauchy problem for the quadratic cascade equation. 
In Section 3 (the Appendix 1) we give the exposition of ultrametric analysis. 
In Section 4 (the Appendix 2) we put the proof of the crucial Lemma [TJ 
In Section 5 (the Appendix 3) we discuss the cascade models of turbulence. 



2 The quadratic cascade equation 

Consider the quadratic cascade equation 
d 



df y(x,t) + J J F(sup(x, £, v)W& *) Mv, t) - v(x, t)) dv{t)dv{ V )+ 

+ J G(svLp(x,t))(v(x,t)-v(Z,t))dv(t;) = 0. (1) 

Here the complex valued functions F and G on the tree T(X) of balls in the (locally compact) 
ultrametric space X with the Borel measure v should satisfy the condition of convergence 

J2\F(J)MJ)-v(JJ))<™ (2) 
j>i 

(analogously for G), where u(J, I) is the measure of the maximal subball in J which contains the 
ball I. 

The third term in the above equation contains the ultrametric pseudodifferential operator, see 
the Appendix 1. The second term is the integral operator with quadratic nonlinearity. 
We have the following important lemma. 

Lemma 1 Let <ft, ip be ultrametric wavelets. Denote J (correspondingly I ) the minimal ball 
which contains the support of ip (correspondingly of <p). 

Let the integration kernel F in (0) satishes condition (0). Then the integral 



J[0,^](x) = J J F(sup(x,e^))0(e)(V(77)-^(x))dKO^) ( 3 ) 
converges and has the form: 

T[(j)^](x) = ipix^x^u, 

where the coefficient $/j can be non-zero only for the case when J < I (the ball J is a strict 
subset of the ball I). In this case 

<b IJ = v\l,J)F{I)-v\j)F{J)- (v 2 (L)-AL,J))F(L). (4) 

L:J<L<I 
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For the proof of the above lemma see Appendix 2. 
Remark If we fix in ((3]) the order of integration in the following way 



F(sup(x, f , 77))0(O mv) ~ 1>(x)) dv(£) 



di/(rj) 



then the integral will converge for any <ft, ^ from D pT) without the additional condition (j2J). 

Example 1 Using lemma [T] one can easily construct the particular solution of the equation 
(j2J) in the form 

v (x,t) = e~ mt ^(x) 

where ip is a wavelet (with the support in the ball I) and 77/ is the corresponding eigenvalue for 
the pseudodifferential operator in the third term of ((2J): 

^ = G{I)v{I) +^G(J)(KJ) - i/(J,I)). 

The introduced function will be a solution of (T3J) because the integral X[ip, xjj], where ip is a wavelet, 
is identically zero. 

Example 2 Let us find a solution of (j2j) in the form 

u(x, t) = ui(t)-0(a;) + v 2 (t)(f)(x), 

where ip an d are wavelets with the supports in the balls J and / correspondingly. 

If / and J have zero intersection or / = J then I[4>, ip] = T[ip, <fi] = and the corresponding 
solution of ([2]) takes the form 

v(x,t) = «i(0)e-^ J V(«) + w 2 (0)e~ ??/ V(x). 

Let the balls / and J be comparable and different, for instance let J < I. Then the corre- 
sponding solution of fl2]) takes the form 

v(x,t) = i; 1 (0)e^ r ' J * + ^ ) *" 1 ' 2(0) ^ 1 ( e ' r ' 7t - 1 )v(a;) + v 2 {Q)e^' t (j){x). 
Here $/j is given by (jl]) and we assume T]i ^ 0. 

We have constructed the nontrivial solution of the nonlinear equation ([2]). This solution 
contains a double exponent (exponent of exponent) of time. 

Let us consider Cauchy problem for equation <^ with v = v(x,t) G D (X) <S> C' 1 ([0, oo)). The 
space Dq(X) £§> C 1 ([0, oo)) is the inductive limit of spaces D Q (S) <S> C 1 ([0, oo)). This means that 
any function v G D (X) ® C 1 ([0, oo)) for any t belongs to some D(S) (where S does not depend 
on t) and v(x,t) is continuously differentiable with respect to t for any fixed x. 

We say that v G D (X) ® C 1 ([0, oo)) is the solution of Cauchy problem for (j2J) with the initial 
condition t> G £>o(X), if 

v(x, 0) = t'o(^) 

and v satisfies ([2]) for t > 0. 

The following theorem describes solutions of Cauchy problems for ([2j) with initial conditions 
mD (X). 
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Theorem 2 Cauchy problem with the initial condition vq G Dq(X) for the integro-differential 
equation (Q3) which satisfies the condition (TJP possesses the unique solution v G Dq(X)<S)C 1 ([0, oo)). 



Proof Any function v G D (X) (g> C 1 ([0, oo)) has the form of the finite linear combination of 
wavelets 

v(x,t) = 5^«/j(t)^/j(a;), 

where Vjj are in C 1 ([0, oo)) and I belong to some S\S m i n , S C T(X). Substituting this into (j2J) 

we get 

r d 
dt 



j( x ) 



0. 



— vij(t) +n I v Ij (t) +v Ij (t) ^ VJj'i^Ji'i^Ji 

Jj':J>I 

Since wavelets are linearly independent and wavelets of larger scale are constants on the supports 
of wavelets of smaller scale, the above equation is equivalent to the system of ordinary differential 
equations 



d , . 



-vjj(t) Vi+ ^2 v J3'i t )' l l J Jv{ x )^Ji ■ (5) 

Jj':J>I 

Here the summation runs over the increasing sequence of balls which are larger than I. This 
system is nonlinear (quadratic). Note that, since x G I and the wavelet ipjj for J > I is constant 
on J, the coefficient ipji 1 ( x ) does not depend on the choice of x G I. 

Let us describe the recurrent procedure of construction of the solution for the system (j5J). The 
initial condition for (121) as a function in D (X) has the expansion over wavelets 



v[x 



0) = $>/,(0>M:r) 



(6) 



Let as choose a maximal / for which the above initial condition vjj(0) is not equal to zero (such I 
can be non unique). Since for a maximal / the corresponding equation in (jS]) is linear we obtain 
the corresponding exponential solution vij(t) of the Cauchy problem: 

««(*)= «/i(0)e-"*. 

Then for maximal subballs /' < I we substitute the obtained solution in the corresponding 
equations in (jSJ) and obtain the equations for vpf. These equations will be linear and will depend 
on the computed at the previous step function vij(t). 

Then we iterate this procedure and get solutions for all pairs 



Vlj (t) = ^(Oje^^fi^w"^'' 1 '^''' 1 ' 



(7) 



It is easy to see that if vjj(0) = then Vjj(t) = 0. Therefore, since the initial condition is a finite 
linear combination of wavelets, the described recurrent procedure will give all non zero vjj(f) 
in a finite number of steps. 

The uniqueness of the solution follows from the observation that any equation in (jSJ) is linear 
as equation for vij(t) and depends only on functions vjj>(t) related to larger scales. This finishes 
the proof of the theorem. □ 
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Remark The formula ([7J gives the explicit solution of Cauchy problem for equation (J2]). 
Therefore this equation is exactly solvable, and the solution can be described with the help of the 
recurrent hierarchical procedure. 

Remark If the initial condition vq for ([2]) belongs to some (finite dimensional) space Dq(S), 
S C T(X), then for any t the solution v(x, t) of the corresponding Cauchy problem for ([2]) will 
also belong to D (S). This observation is the example of the mentioned in [S] general phenomenon 
of existence of localized solutions for some integral equations in ultrametric analysis. 

Remark The equation (T2J) for v in D (X) ®C 1 ([0, oo)) is equivalent to the system of ordinary 
differential equations (jSJ). The system (jSJ) gives the example of the so called cascade model. Models 
of this kind are used for description of turbulence, see the Appendix 3. The relation of nonlinear 
ultrametric integral equations and cascade models of turbulence was mentioned for the first time 
in [8]. For discussion of cascade models of turbulence see [TO] . 

The results of the present paper show that the cascade models are related to integrable non- 
linear ultrametric integral equations and ultrametric wavelet analysis. 

3 Appendix 1: Ultrametric analysis 

In this Section we discuss some results on ultrametric analysis, which can be found in [5], [6], [TJ. 

Definition 3 An ultrametric space is a set with the ultrametric d(x, y) (where d(x, y) is 
called the distance between x and y), i.e. a function of two variables, satisfying the properties of 
positivity and non degeneracy 

d(x,y)>0, d(x,y) = ==► x = y; 

symmetricity 

d(x,y) = d(y,x); 

and the strong triangle inequality 

d(x, y) < max(d(x, z), d(y, z)), Wx, y, z. 

We say that an ultrametric space X is regular, if this space satisfies the following properties: 

1) The set of all the balls of nonzero diameter in X is finite or countable; 

2) For any decreasing (infinite) sequence of balls {I^ k '}, I^ k ' D I^ k+1 \ the diameters of the 
balls tend to zero; 

3) Any ball of non-zero diameter is a finite union of maximal subballs. 

Ultrametric spaces are dual to directed trees. Below we describe some part of the duality 
construction. 

For a regular ultrametric space X consider the set T(X), which contains all the balls in X of 
nonzero diameters, and the balls of zero diameter which are maximal subbals in balls of nonzero 
diameters. This set possesses a natural structure of a directed tree. Two vertices / and J in T(X) 
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are connected by an edge if the corresponding balls are ordered by inclusion, say / D J (i.e. one 
of the balls contain the other), and there are no intermediate balls between I and J. 

The partial order in T(X) is defined by inclusion of balls, this partial order is a direction. We 
recall that a partially ordered set is a directed set (and a partial order is a direction), if for any 
pair of elements there exists the unique supremum with respect to the partial order. 

On the directed tree T(X) we have the natural increasing positive function which puts into 
correspondence to any vertex the diameter of the corresponding ball. 

Assume now we have a directed tree T with the positive increasing function F on this tree. 
Then we define the ultrametric on the set of vertices of the tree as follows: d(I, J) = F(sup(J, J)) 
where sup(J, J) is the supremum of vertices /, J with respect to the direction. 

Then we take completion of the set of vertices with respect to the defined ultrametric and 
eliminate from the completion all the inner points of the tree (a vertex of the tree is inner if it 
does not belong to the border of the tree). We denote the obtained space X(T), this space is 
ultrametric. 

An ultrametric pseudodifferential operator is defined in the following way. Consider a a- 
additive Borel measure v with countable or finite basis on a regular ultrametric space X. Consider 
the pseudodifferential operator 



Here T(I) is some complex valued function on the tree T(X). The supremum sup(x, y) — I of 
the points x, y G X is the minimal ball / in X, containing both points. 

Let us build a basis in the space L 2 (X, v) of complex valued functions on a regular ultrametric 
space X which are quadratically integrable with respect to the measure v. We call this basis the 
basis of ultrametric wavelets. 

Denote by V(I) the space of functions on X, generated by characteristic functions of the 
maximal subballs in the ball I of nonzero diameter. Correspondingly, Vo(/) is the subspace of 
codimension 1 in V(I) of functions with zero mean with respect to the measure v. The spaces 
V (I) for different / are orthogonal. We denote pj the number of maximal subballs in the ball I. 
The dimension of V (I) as the Euclidean space with the scalar product as in L 2 (X, v) will be less 
or equal to pi — 1. 

We introduce in the space Vq(I) some orthonormal basis {ipij}- If the measures of all maximal 
subballs in I are positive, the index j can take values 1, . . . ,pi — 1. The next theorem shows how 
to construct the orthonormal basis in L 2 (X, v), taking the union of bases {ipij} in the spaces V (I) 
over all non minimal I G T(X) (equivalently, over all balls / of non-zero diameters). 

Theorem 4 1) Assume that the measure v(X) of the regular ultrametric space X is infinite. 
Then the set of functions {ipij}, where I runs over all non minimal vertices of the tree T(X) is 
an orthonormal basis in L 2 (X, v) . 

2) Let the measure v{X) of the regular ultrametric space X is finite and is equal to v(X) = A. 
Then the set of functions {ipij, A~2 } ; where I runs over all non minimal vertices of the tree T(X) 
is an orthonormal basis in L 2 (X, v) . 

The basis introduced in the present theorem will be called the basis of ultrametric wavelets. 
The next theorem shows that the basis of ultrametric wavelets is the basis of eigenvectors for 
ultrametric pseudodifferential operators. 
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Theorem 5 Let the following series converge absolutely: 

£t(J)(i/(J)-i/(J,J2)) (8) 

J>R 

for some ball R. 

Then the ultrametric pseudodifferential operator 

Tf(x) = J T(sup{x,y)){f{x) - f{y))du{y) 

has a dense domain in L 2 (X, v) and ultrametric wavelets from Theorem^ are eigenfunctions ofT: 

T^/j(x) = X I i[) Ij (x) (9) 

with the eigenvalues: 

X T = T(I)u(I) + J2 T(J)HJ) ~ »{J, I))- (10) 
j>i 

Here (J, I) is the maximal subball in J which contains I. 

A function / on an (ultrametric) space X is called locally constant, if for any arbitrary point 
x G X there exists a positive number r (depending on x), such that the function / is constant on 
the ball with the center in x and the radius r: 

f(x) = f(y), Vy: d(x,y) <r. 

In particular, the characteristic function xi of a ball I is locally constant. 

The space of test functions D(X) on the ultrametric space X is defined as the space of lo- 
cally constant functions with compact support. Any test function in D(X) is a (finite) linear 
combination of characteristic functions of balls. 

The space D(X) is the inductive limit of the finite dimensional spaces D(S), where S is a 
regular subtree in T(X) and D(S) is the linear span of characteristic functions of balls from S. 
Here the finite subtree S C T(X) is called regular, if for any increasing edge IJ, I, J E S, I < J 
the subtree S contains all edges KJ, K < J. 

We denote Vq(S) C D(S) the subspace of mean zero functions in D(S). We denote D (S) the 
factor space Vo(S)/V^-(S) over the subspace Vq~(S) C Vo(S) of functions which can be not equal 
to zero only on some sets of zero measure v. The space Dq(S) can be identified with the linear 
span of ultrametric wavelets ipij, I £ <S\«S m i n , where iS min is the set of the minimal elements in S. 

We denote D Q (X) the inductive limit of the finite dimensional spaces D Q (S). The space D (X) 
is isomorphic to the linear span of ultrametric wavelets. 



4 Appendix 2: Proof of Lemma H 

Consider the integral 

J[0,^](x)= [ [ F( S u P (x,£, V ))^)(i;( V )-i;(x))dv(Z)dv(v) 
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for the case when and ip are ultrametric wavelets. The condition ([2]) guarantees the convergence 
of the integral. Let us denote I and J the minimal balls which contain the supports of and ip 
corresp ondingly. 

Consider the following cases: 

1) The balls J and / have empty intersection. Then the integral over £ gives zero because is 
a mean zero function and the integration kernel F is constant on /. 

2) Let I be the strict subset of J. Then the integral over £ can be nonzero only for x,rj G /. 
In this case ip(rj) and ip(x) cancel and thus the integral is equal to zero. 

3) Let I = J. Then the integral over £ can be nonzero only for x, rj G /, and, moreover, 
sup(x, rj) < I. In this case ip(rj) and ip{ x ) cancel and thus the integral is equal to zero. 

4) Let J be the strict subset of I. Consider the following cases. 

a) If x ^ J then ip(x) = and the integration kernel is constant on J. Therefore the integral 
over rj is equal to zero. 

b) Let x G J < I. The integral takes the form 



T[<j>,1>](x) = / / F(sup(a;,f,77))0(O (Vfo) - Mv)M0 



F(sup(x, f , ??))0(£) (^(77) - ip(x)) dv{rj)dv{i) 



££l J rj:sup(x,r))<I 



£:sup(x,£)=I J §:sup(:r,{;)</J J ri:sup(x,ri)<I 

= -cf>{x)v(I, J)F(I) 



F(sup(x, f , 17)) 0(0 (ip(rj) - if>(x)) dv(7])dv(£) 
(ip(rj) - 4>(x)) du(rj)+ 



+0(x) 



+ 

£:sup(x,£)<J J £:J<sup(z,§)</ 



ri:sup(x,ri)<I 

F(sup(a;, f , 7/)) (ip(r]) - ip(x)) dv{rj)dv{£) 

r]:sup(x,ri)<I 



^j(x)(j)(x)u 2 (I,J)F(I) + (f)(x)u(J) / F(sup(x,rj))(ij(ri) -ip(x))dv(rj)- 

' rj: J<sup(x,rj)<I 



-ip(x)(p(x) / / F(sup(x,£,rj))dis{rj)dis{£) 

J §: J<sup(x,l;)<I J r):sup(x,r))<I 



ijj{x)4>{x) 



Xl,J)F{I)-v\J)F{J)-v{J) F{L){v{L)-v{L,J))- 



L:.J<L<I 



(v(L)-is(L,J))lF(L)is(J) + F(sup(L,K)) ^ (u(K) - u(K, J)) 

\ K:J<K<I , 

(I,J)F(I)-S(J)F(J)-2v(J) F{L){v{L)-u{L,J))- 

L:J<L<I 

J2 (y(L)-u(L,J)) J2 HK)-u(K,J))F(sup(L,K)) 



L:J<L<I 

= ip(x)<p(x) 



L:J<L<I 



K:J<K<I 
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Since 



(v(L)-u{L,J)) HK)-u(K,J))F{Bap{L } K)) = 




= £ (v(L)-v(L,J)) 2 F(L) + 2 {v{L)-v{L,J))F{L){v{L,J)-v{J)), 



L:J<L<I L:J<L<I 

we get for the integral 



T[M\{x) = ^{x)(j>{x) v 2 {I,J)F{I) 



v\J)F{J) 



(v(L)-v(L,J)) 2 F(L)-2 (v(L)-v(L,J))F(L)v(L,J) 




1>{x)<Kx) v 2 {I,J)F{I)-v 2 {J)F{J)- (v 2 (L)-v 2 (L,J))F(L) 



L:J<L<I 



This finishes the proof of the lemma. □ 

5 Appendix 3: Hierarchical models of turbulence 

In the present Section we discuss, following the book [ID] , the hierarchical (or cascade, or shell) 
models of turbulence. The construction [TUJ of cascade models is as follows. 

1) The hierarchical basis in L 2 (M 3 ), similar to bases of wavelets, is constructed. The hierarchical 
basis is built using translations and dilations of some finite set of mean zero functions. 

2) The hierarchical basis is applied to the investigation of the Navier-Stokes equation 



Here v is velocity, p is density, p is pressure, fi is viscosity, and f is external force. 

This means that we substitute in the above equation the expansion of velocity over the hierar- 
chical basis, and investigate the obtained system of non-linear ordinary differential equations for 
the coefficients of the expansion: 



3) We make the following important observation: in the obtained at the previous step system of 
ordinary differential equations majority of the coefficients will be small, and, moreover, will decay 
with the distance from the main diagonal, i.e. A^u decays with \N — I\ —>■ oo or | N — J\ — » oo, 
Bnk decays with \N — K\ — > oo. 

We define the cascade system of ordinary differential equations as follows: we put in ( fl2l) the 
coefficients Anij, Bnk, corresponding to sufficiently large \N — I\, \N — J\, \N — K\ equal to 
zero. We obtain the system of equations with contributions corresponding to finite number of 
main diagonals of the tensors A^u, B^k- In particular, any equation in the obtained system will 




(11) 



dv N (t) 

dt 



A mJ Vi{t)vj{t) + B NKV K (t) + C N . 



(12) 



IJ K 
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contain finite number of contributions. In this sense the cascade model is constructed by cut-off 
of the system (fT2l . 

The obtained in this way cascade systems of differential equations were studied numerically 
[TU] . It was found that the behavior of some cascade models is in good agreement with the 
experimental results. 

The hierarchical structure of the cascade models is well suited for the description of the Richard- 
son cascade and self-similarity of turbulent flows [10] . The indices of the hierarchical basis in (fT2l) 
correspond to eddies with different scales and positions. Tensor A^u describes interaction of 
eddies, tensor B^k describes dissipation of energy due to viscosity. 

Let us discuss the relation of cascade models and ultrametric analysis. The natural conjecture 
is that the cascade system of ordinary differential equations is equivalent in some sense to the single 
integral equation with integration over the ultrametric argument. The first result of this kind was 
obtained in [8], where the p-adic integral equation which describes the Richardson cascade was 
constructed. 

In the present paper we introduce the approach related to ultrametric wavelet analysis. We 
built the integro-differential equation (j2J) with the following property: the system (j5J) obtained by 
expansion of v(x,t) over the basis of ultrametric wavelets is an example of the cascade system of 
differential equations. In this approach the eddies correspond to ultrametric wavelets, which itself 
are in correspondence with the elements of the hierarchical basis of [10] . The second (quadratic) 
term in (j2J) describes interaction of eddies, the third term in (J2J) describes dissipation. The cut- 
off procedure used in the construction of the cascade model in our approach takes the form of 
approximation of the Navier-Stokes equation (TTTj) . in the regime of developed turbulence, by the 
cascade equation The approximation here is understood as approximation of the system (|T2l 
by the system ([5]). 

Actually in the present paper we have described by means of ultrametric methods the simplest 
cascade model. We found that this simplest cascade model is integrable. This shows the relation 
between cascade models and theory of integrable nonlinear integral ultrametric equations. 
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